The equilibrium state of a trapped two-dimensional Bose gas 
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We study experiment ahy and numericahy the equihbrium density profiles of a trapped two- 
dimensional ^^Rb Bose gas, and investigate the equation of state of the homogeneous system using 
the local density approximation. We find a clear discrepancy between in-situ measurements and 
Quantum Monte Carlo simulations, which we attribute to a non-linear variation of the optical den- 
sity of the atomic cloud with its spatial density. However, good agreement between experiment and 
theory is recovered for the density profiles measured after time-of- flight, taking advantage of their 
self-similarity in a two-dimensional expansion. 
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Low-dimensional atomic gases provide stringent tests 
of the many-body description of quantum matter, be- 
cause thermal or quantum fluctuations play a more im- 
portant role than for three-dimensional (3d) fluids [Ij. 
These systems are prepared by freezing one or two mo- 
tional degrees of freedom [2j. For 2d Bose gases, recent 
experiments [3 , 4J and corresponding numerical analyses 
O [6| gave evidence for a Berezinskii-Kosterlitz-Thouless 
transition, with a quasi-long range order of the phase of 
the gas below a critical temperature. 

A remarkable feature of the uniform 2d Bose gas is the 
scale invariance of its equation of state. For a large do- 
main of parameters, the phase space density D = nX^ 
is not an independent function of the chemical poten- 
tial /i and the temperature T, but depends only on the 
ratio a = ji/k^T. Here n is the 2d spatial density 
and A = (27rfi^/m/cBT)^/^ is the thermal wavelength. 
Scale invariance originates from the fact that the in- 
teraction strength in 2d is determined by a dimension- 
less number that is approximately energy-independent, 
g = a^dl where asd is the scattering length char- 
acterizing low energy interactions in 3d, and iz is the 
thickness of the gas along the frozen direction z [for har- 
monic confinement with frequency cj^, — {h/muJzY^'^] 
Since interactions provide no energy scale, dimen- 
sional analysis implies that D can be cast in the form 
D = F{a^g). This holds when £z ^ <^3d, and applies for 
all experiments so far. 

In this Letter we present results from combined experi- 
mental and numerical studies to test this scaling property 
of the 2d Bose gas. The experiments are performed with 
^^Rb atoms and the results are compared with Quantum 
Monte Carlo (QMC) simulations. Experimental density 
distributions are inferred from the absorption of a reso- 
nant probe light beam. We investigate both the in-situ 
distribution of the gas and the one obtained after a time- 
of-flight (TOF) in the xy plane. The numerical results 



confirm the scale invariance and are in excellent agree- 
ment with the prediction based on the equation of state 
for the uniform 2d gas [8j and on the local density ap- 
proximation. The measured in-situ distributions clearly 
differ from the numerical predictions if we assume a lin- 
ear relation between the optical density of the cloud and 
its spatial density n. However we point out that the usual 
single-scattering approximation for the probe beam pho- 
tons, which is at the basis of this linear relation, is insuf- 
ficient in our situation. Indeed the inter-particle distance 
in the center of the trap is comparable to /c^ ^, where /cl 
is the wave vector of the probe light. After long TOF 
durations the densities are considerably lower, and the 
single scattering approximation holds. We then recover 
good agreement between theory and experiment, using a 
second (dynamical) scale invariance of the 2d Bose gas 
confined in a harmonic potential [9j. 

In our experiment, we first prepare a Bose-condensed 
gas of - 3 X 10^ ^'^Rb atoms in the F = 2, mp = 2 
hyperfine state of the ground level 5^Si/2- The gas is 
confined in a magnetic Time-averaged Orbiting Potential 
trap [To] at a temperature of 160 nK, obtained through 
evaporative cooling with a radio- frequency (rf) field. We 
then add a dipole potential providing strong confinement 
in the vertical direction. This potential is generated by 
a laser beam at a wavelength of 532 nm and a power 
of 1.3 W. The beam passes a holographic plate that im- 
prints a phase of tt on its upper half, and is focussed 
onto the atoms [llj. There it creates (together with the 
magnetic trap) the potential illustrated in Fig. [l^. The 
vertical and horizontal waists of the beam in the absence 
of the holographic plate are 5.0 |LLm and 140 iitm, respec- 
tively. The trapping frequencies in the combined poten- 
tial, measured by exciting the dipole oscillation of the 
atom cloud, are ujz/27r = 3.6 (3) kHz along the vertical 
axis, leading to ^ = 0.146(6), and uJxI^h = 21.0 (5) Hz 
and ujy /27r = 18.8 (5) Hz in the horizontal plane. For 
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FIG. 1: (Color online) (a) Potential V along the vertical di- 
rection z created by the magnetic trap and the laser beam, 
(b-c) Side view of the cloud before (b) and after (c) depump- 
ing atoms in the side wells. The horizontal stripes are due to 
diffraction, (d) Top view in-situ image yielding fit parameters 
T = 132 nK and a = 0.29 (for ^ = 0.25). 

convenience we define uj = (cja^cj^)^/^ = 27r x 19.9 Hz. 

The dipole potential is ramped up in 1.5 s and it 
cuts out a plane of atoms in the center of the original 
cloud, while the remaining atoms reside in the side wells 
(Fig. We depump these unwanted atoms to a non- 
detected state (F = 1) with a 35 ms laser pulse tuned to 
the F = 2 ^ F' = 2 transition of the D2 line, masking 
the central plane of atoms (Fig. [l]:). The atoms are held 
in the combined trap during 5 s in the presence of the rf 
which controls the final temperature. 

We probe the spatial density distribution in the xy 
plane using absorption imaging of a resonant probe laser 
beam propagating along the vertical axis. The mea- 
surement provides a map of the optical density (OD) 
of the atomic cloud, defined as the natural logarithm of 
the ratio between incident and transmitted light intensi- 
ties. Neglecting multiple scattering of photons, one has 
OD(r) = crn(r), where a is the absorption cross-section. 
Absorption imaging can be performed in the presence of 
the trap magnetic field (in-situ measurement), or after 
a TOF expansion in the xy plane. A typical in-situ im- 
age is shown in Fig. [l]i. From each image we generate 
a radial density profile OD(r) by averaging over the az- 
imuthal angle while accounting for the residual ellipticity 
in the xy plane. 

We determine ji and T by fitting the profile OD(r) 
with the mean-field Hartree-Fock (MFHF) prediction, 
which reads for a strictly 2d gas (see e.g. [12j): D = 
— In (1 — e'^-gD/i^^^ Since our highest temperatures (~ 
150 nK) are comparable to hjOz/k^ = 170 nK, the gas 
is only quasi- 2d [13-r.-15j and we take into account cor- 
rections to D{a) due to the residual excitation of the 
atom motion in the z direction ^3 We restrict 
the fitting domain to simultaneously fulfill two condi- 
tions: (i) the phase space density in the axial ground 
state must be lower than 2.5, so that beyond-mean- field 
corrections are negligible [13, 16j; (ii) the optical density 
of the cloud must be lower than 0.2 to exclude distortions 
due to multiple scattering of probe photons (see below). 
Because the spatial density enters non-linearly in the re- 
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FIG. 2: (Color online) (a) Dots: Measured in-situ profile 
OD(r). For ^ = 0.25 the fit with MFHF theory yields T = 
126 (6) nK, a = 0.34(9), where uncertainties represent stan- 
dard deviations obtained by fitting individual images. Con- 
tinuous line: corresponding QMC simulation with N = 73900 
atoms (inset: same data in log plot). Upper (dash-dotted) 
and lower (dashed) lines: QMC results obtained assuming 
^ = 0.21 [fit parameters (T (nK), a, Ar) = (130, 0.39, 96300)] 
and ^ = 0.29 (122,0.29,57900), respectively, (b) Set of mea- 
sured density profiles (dots) and corresponding QMC results 
(lines) for rf evaporation parameters. From bottom to top: 
(T (nK), a, Ar)=(87, 0.49, 54100) [black], (109,0.39,63800) 
[red], (142, 0.28,78400) [blue], (153, 0.23,79900) [magenta]. 
Each experimental profile in (a) and (b) is an average of 9 
images. 



lation between D and a, the fitting procedure can also 
provide a value of the detection efficiency ^, defined as 
the ratio between the actual absorption cross-section and 
the ideal one expected for monochromatic probe light in 
the absence of stray fields. For all images corresponding 
to a given TOF duration we extract a single value of ^ 
from the fit. For in-situ images we find the small value 
^ = 0.25 (4), which accounts for the strong reduction of 
the absorption cross-section due to the magnetic field of 
the trap. With this correction factor, an optical density 
of 0.1 corresponds to a density n c:^ 3.0 \xm~^. 

Our path-integral QMC simulations are performed in 
the canonical ensemble in 3d continuum space, for the 
same geometry as the experiments. They take into ac- 
count residual excitations along the strongly confined di- 
rection z [5] [17]. Pair interactions are described by a 3d 
pseudo-potential [l8j. All thermodynamic properties of 
the gas are obtained to high precision and without sys- 
tematic errors. The chemical potential associated to a 
given atomic distribution is obtained from a fit of the 
wings of the distribution with the MFHF prediction, as 
for experimental data. 

We now compare our experimental and numerical re- 
sults, first confronting the measured in-situ optical den- 
sities with theoretical profiles calculated for the same fi 
and T. As illustrated in Fig. [2] (a-b), the wings of the 
calculated and measured profiles nearly coincide. How- 
ever there is a clear discrepancy in the central part of 
the density distributions. Whereas the central OD in 
the four coldest experimental distributions is ~ 1.0, the 
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FIG. 3: (Color online) Combination of the experimental and 
theoretical results of Fig. [2] Error bars indicate the standard 
deviation of the data, (a) Measured OD as a function of calcu- 
lated OD, averaged over the data shown in Fig. [2] The dashed 
line with a slope 1 is a guide to the eye. (b) Continuous lines: 
QMC results for D as a function of aiocai = a — muJ^r^ jk^T 
for the data shown in Fig. |2] (same color code). Black dashed 
line: prediction of [8 for the uniform case. Dots: Measured 
averaged over all experimental data shown in Fig. |2] 

QMC simulations systematically predict a central OD 
~ 1.8, i.e. a density of 55 jitm"^. A global com- 
parison between predicted and measured OD is shown 
in Fig. [3^, where we performed an average over the 5 
profiles of Fig. [2j 

We now discuss several possible causes for this discrep- 
ancy. A first possible source of error is the uncertainty 
on the detectivity factor ^. To estimate its influence, we 
have reprocessed the measured proflle shown in Fig. [2^, 
by choosing the lower (^ = 0.21) and upper (^ = 0.29) 
values of the uncertainty interval for ^. The QMC re- 
sults for the modifled flt parameters are shown with dash- 
dotted and dashed lines in Fig. |2^. Clearly, the uncer- 
tainty on ^ does not account for the observed deviation. 
Other 'technical' causes for this discrepancy could be the 
imperfect resolution of the imaging system and/or the 
atomic motion during the imaging pulse. However, nei- 
ther of them can account for the difference between pre- 
dicted and measured density proflles [19j. The most prob- 
able cause of this discrepancy is the reduced absorption 
cross- section for large 2d atomic densities, due to mul- 
tiple scattering of the photons of the probe laser beam. 
Although our optical densities (< 1) do not exceed usual 
values for absorption imaging, they correspond in this 2d 
geometry to a short mean distance d between scatterers. 
For the densest clouds we flnd that k\,d is on the order 
of 1 (/cl = 8 X 10^ m~^), which can signiflcantly modify 
the photon scattering rate [20l |2TJ . 

The trapped system of interest here can be related 
to a uniform gas using the local density approximation 
(LDA). Within LDA, the phase space density D(r) in 
a trapping potential V{r) is given by D(r) = F(/i — 
F(r),T) where D = F{ii^T) is the equation of state of 
the uniform system. For the quasi-2d Bose gas, the va- 
lidity of LDA has been accurately checked with QMC 
simulations in [17J. To test the scale invariance of the 
equation of state, we have plotted in Fig. the phase 
space density I) as a function of the ratio aiocai = ocir) 



between /i — ma;^r^/2 and k^T^ using the same data as 
in Fig. [2J The calculated functions D{a(r)) (continu- 
ous lines) nicely superpose onto each other, conflrming 
the scale invariance in this low temperature region where 
the excitation of the z motion does not play an impor- 
tant role. Note that the central phase space densities 
notably exceed the critical value for superfluid transition 
{Dc ^ 8.0 for g = 0.146 [22j), signaling the presence of a 
signiflcant superfluid component. The QMC results con- 
flrm the prediction of [8j for the uniform 2d gas, which 
was obtained using a classical fleld Monte Carlo method 
(dashed line in Fig. The small corrections to the 

results of [8] are due to two factors: (i) the presence of 
residual excitations in the z-axis, and (ii) the flnite value 
of the interaction parameter g [17 \. Obviously the devia- 
tion between experimental and numerical data that was 
appearing in the density proflles of Fig. |2] also shows up 
in the plot of D{a{r)) in Fig. |3]3, and the experimental 
values of D (dots) lie systematically below the predicted 
ones in the high phase space region. 

A simple way to circumvent the problem of imag- 
ing high density regions is to take advantage of another 
known scale invariance of the 2d Bose gas. This invari- 
ance manifests itself in a 2d ballistic expansion after the 
conflnement in the xy plane has been suddenly released. 
It follows from a hidden SO (2,1) symmetry of the 2d 
Bose gas with contact interactions V{r) = {h^g /m)5{r) 
when it is conflned in an isotropic harmonic potential of 
frequency uo [9j. Starting from an arbitrary initial equi- 
librium proflle neq(r), the density proflle after a TOF 
duration t is obtained by a scaling transform 

n(r, t) = T]^ ne^{r]tr) , = (1 + ujH^)-^^^ . (1) 

This relation was predicted within the Bogoliubov ap- 
proximation in [23j and it holds exactly for interaction 
potentials that satisfy V{Xr) = V{r)/X'^ [9]. The 2d con- 
tact potential belongs to this class of functions, although 
strictly speaking it needs to be regularized at short dis- 
tances to avoid ultraviolet divergences. 

Experimentally, we initiate the 2d expansion by switch- 
ing off the magnetic trap while keeping the optical po- 
tential constant. The atom cloud thus expands in the xy 
plane for an adjustable duration t after which we take 
an image of the cloud. We explore TOF durations up 
to t = 14 ms for which the central density is reduced by 
a factor 7^^~^ = 4, so that artifacts due to multiple scat- 
tering of probe photons should be strongly reduced (see 
Fig. [3^). We show in Fig. [4^ a succession of density pro- 
flles recorded for TOF durations varying from to 14 ms. 
Each proflle has been rescaled to the initial in-situ distri- 
bution according to the law ([T]) , so that ideally all proflles 
should be superimposed. In practice this superposition is 
poor for short TOF durations and becomes better as the 
clouds expand. All scaled proflles obtained for t > 10 ms 
coincide within their noise, showing that the detectivity is 
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FIG. 4: (Color online) (a) Optical density profiles obtained 
for a TOF duration t = (black), 3 (red), 6 (green), 10 (blue), 
12 (cyan), 14 (magenta) ms and rescaled to their in-situ value 
according to (IT]) . (b) Squares: Optical density profile obtained 
by averaging the results of (a) for 10 < t < 14 ms, yielding 
fit parameters T = 94 nK, a = 0.36 for ^ = 0.63. Lines: 
QMC results for the same fit parameters (continuous, N = 
42000), and for those deduced assuming ^ = 0.47 [dashed red, 
(T (nK), a, N)= (104, 0.39, 57600)] and ^ = 0.79 [dash-dotted 
blue, (87,0.33, 32100)]. 

constant. We find ^ = 0.63 (16), consistent with the mea- 
sured absorption line- width of the probe laser. Moreover, 
the agreement between these scaled profiles and the cal- 
culated ones is indeed very satisfactory (see Fig. [4]3)[24j. 
Note that the scale invariance accessible with such 2d ex- 
pansions is notably different from the properties revealed 
in a 3d expansion of an initially quasi-2d gas, for which 
all trapping potentials are switched off simultaneously. 
In that case interactions are negligible during the TOF, 
due to the fast expansion along the initially strongly con- 
fined direction. This gives access to the quasi-coherent 
core which is present at the center of the trap below the 
Berezinskii-Kosterlitz-Thouless transition |3j [H [5^ . 

In conclusion we have studied the equilibrium density 
profile of a trapped quasi-2d Bose gas. The discrepancy 
between in-situ absorption images and numerical calcula- 
tions suggests that multiple scattering of the probe pho- 
tons reduces the absorption cross-section in high-density 
regions [21j. To reveal the undistorted density distri- 
bution of the gas, we have taken advantage of its self- 
similarity in a 2d expansion. The profiles measured af- 
ter this expansion are in good agreement with the QMC 
predictions. Natural extensions of this work are the mea- 
surement of other thermodynamical quantities from in- 
situ images using the procedure proposed in [26], and 
the study of the regime where the scale invariance breaks 
down and energy-dependent corrections to g become im- 
portant [T6[l27H3Q]. 
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